
subjectList = {...
    'A080512A'...
    'K091412S'...
    'M101012R'...
    'N072912R'...
    'R091412V'...
    };

numVox = 200;


for subjIndx = 1:length(subjectList);
    
    [allResults, averageBetas] = extractROImask(subjectList{subjIndx}, numVox);
    
    display(averageBetas);
    
    figure(1);
    subplot(length(subjectList),1,subjIndx);
    bar(allResults{3,4});
    title(allResults{1,1});
    ylabel('Average Beta value');
    hold on;
    
    figure(2);
    subplot(length(subjectList),1,subjIndx);
    bar(allResults{2,4});
    title(allResults{1,1});
    ylabel('Average Beta value');
    hold on;
    
    figure(3);
    subplot(length(subjectList),1,subjIndx);
    bar(allResults{4,4});
    title(allResults{1,1});
    ylabel('Average Beta value');
    hold on;
    
end

xlabel('Scan number');
